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1. Introduction 

In these proceedings we follow up on the line of work on 2+1 flavor QCD thermodynamics 
with improved staggered fermions by the HotQCD collaboration. The setup of simulations and 
some preliminary results have been reported earlier, e.g. |Jlj, ||, ||, ^ and the continuum extrapolation 
for chiral T c at the physical quark masses is presented in Full data set includes several lines 
of constant physics down to the light quark mass m; = m s /20 for asqtad and m/ = m s /40 for 
HISQ/tree 1 . The lattice spacings cover the range of temperatures T = 130 — 440 MeV with N x = 6, 
8 and 12 for HISQ/tree and T = 148 - 304 MeV with = 8 and 12 for asqtad. 

Performing the continuum limit requires control over cutoff effects. In improved staggered 
discretization schemes the leading 0(a 2 ) errors at low temperature (coarse lattices) originate from 
violations of taste symmetry, that distort the hadron spectrum. For this reason it is important 
to perform simulations on fine enough lattices (large N z in finite-temperature setup) and/or use 
actions with the smallest discretization effects. Analysis of the discretization effects for asqtad and 
HISQ/tree used in this study is presented in Ref. 

2. Chiral transition 

For vanishing light quark masses there is a chiral phase transition which is expected to be 
of second order and in the 0(4) universality class [f7|]. However, universal scaling allows to define 
pseudo-critical temperatures for the chiral transition even for non-zero light quark masses, provided 
they are small enough. For staggered fermions that preserve only a part of the chiral symmetry 
there is a complication: in the chiral limit at finite lattice spacing the relevant universality class is 
0(2) rather than 0(4). Fortunately, in the numerical analysis the differences between 0(2) and 
0(4) universality classes are small so when referring to scaling we will use the term 0(N) scaling. 
Previous studies with the p4 action provided evidence for 0(N) scaling Similar analysis for 

the asqtad and HISQ/tree action establishing if the 0(N) scaling is applicable is performed in [j|] 
and explained below. 

The order parameter for the chiral transition is the chiral condensate 

M^^Spi. (2.1, 

Its temperature and quark mass dependence near the critical temperature can be parametrized by a 
universal scaling function fc and a regular function fM.reg that describes corrections to scaling: 

M b (T,m l ,m s ) = h i / s f G (t/h 1 /P s )+f M ^(T,H), t = 1 ( ^=^Y h = ±H, H=^ (2.2) 

to \ y f u J ho m s 

and r c ° is the critical temperature in the chiral limit. The pseudo-critical temperature can be defined 
as the peak position of the chiral susceptibility 

Xm < l = ^ {W)l (23) 

'The lightest mass m; = m s /40 is used for the analysis of the chiral condensate and susceptibility while all other 
quantities are calculated at m/ = m s /20 for HISQ/tree. 
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Figure 1: An example of a simultaneous fit to the chiral condensate (left) and susceptibility (right) for 
HISQ/tree on N T — 8 lattices. Open symbols indicate the range included in the fit. Dotted black line is an 
extrapolation to the physical light quark mass. 

whose scaling behavior is also described by fc and fM.reg as 

Xm.l T 2 f 1 , 1/5—1 ^ / \ i dfM,reg(T,H)\ 1 Z , , t 

r^ = ^U /z(Z) + dH J' fz®=s\fo(z)-pf {z)], z = - m . (2.4) 

The singular function fc is well studied in spin models and has been parametrized for 0(2) 
and 0(4) groups. For the regular part we consider leading-order (linear) dependence in H and 
quadratic in T: 



( T-T c ° 

jM,reg(T,H) = I a + a\ — h<32 




(2.5) 



Then we are left with 6 parameters to be determined from fitting the data, T®, to, ho, ao, a\ and a^- 
We perform simultaneous fits to M\, and Xm,l f° r the asqtad action on Af T = 8, 12 and the 
HISQ/tree action on = 6, 8 and 12. An example of such a fit for HISQ/tree, N z = 8 is shown in 

FigJ 

Then, performing a combined l/N 2 extrapolation of T c values obtained with the asqtad and 
HISQ/tree action as shown in Fig. ^| we obtain 



r c = (154±8±l)MeV, 



(2.6) 



where the first error is from the fit and the second is the overall error on the lattice scale determina- 
tion. The fits for asqtad and HISQ/tree are constrained to have a common intercept. (See Ref. [||] 
for more details on the fitting procedure and analysis of systematic errors.) To present a combined 
error we add the two errors giving our final value T c = 154 ± 9 MeV. 



3. Deconfinement aspects of the transition 

The deconfinement phenomenon in pure gauge theory is governed by breaking of the Z(N C ) 
symmetry. The order parameter is the renormalized Polyakov loop, obtained from the bare Polyakov 
loop as 



L ren (T) = z(P) N *L bare (P) = z(P) Nt ( ^Tr n U (x ,x) ) , (3.1) 
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Figure 2: Example of a combined fit to the asqtad 
and HISQ/tree data. The fit is constrained to have a 
common intercept at 1 /N% —> 0. 
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Figure 3: The renormalized Polyakov loop in 
pure gauge 51/(2), 5C/(3) and full QCD with asq- 
tad and HISQ/tree. 



where z(j3) = exp(— c(/$)/2). c(j8) is the additive normalization of the static potential chosen such 
that it coincides with the string potential at distance r = 1.5ro with ro being the Sommer scale. 

In QCD Z{N C ) symmetry is explicitly broken by dynamical quarks, therefore there is no obvi- 
ous reason for the Polyakov loop to be sensitive to the singular behavior close to the chiral limit. 
Indeed, the temperature dependence of the Polyakov loop in pure gauge theory and in QCD is quite 
different, as one can see from Fig. 0. Also note, that in this purely gluonic observable there is very 
little sensitivity (through the sea quark loops) to the cut-off effects coming from the fermionic sec- 
tor. While losing the status of the order parameter in QCD, the Polyakov loop is still a good probe 



of screening of static color charges in quark-gluon plasma [10, 11]. 

Other probes of deconfinement are fluctuations and correlations of various charges that can 
signal liberation of degrees of freedom with quantum numbers of quarks and gluons in the high- 
temperature phase. Here we consider quadratic fluctuations and correlations of conserved charges: 



Xi(T) 1 d 2 \nZ(T,^ 



T 2 T^V d(Hi/T) 2 



Z\j(T) _ 1 d 2 hLZ(T,iu,iij) 
*=<,' T 2 T 3 V dijli/T)d(jlj/T) 



(3.2) 

;Ui=M/=0 



Fluctuations are also sensitive to the singular part of the free energy density. The light and strange 
quark number susceptibilities are often considered in connection with deconfinement. These quan- 
tities show a rapid rise in the transition region and for the strange quark number susceptibility the 
behavior is consistent with the Hadron Resonance Gas (HRG) model, as shown in Fig. ||. For the 
light quark number susceptibility, a quantity dominated by pions and therefore very sensitive to 
the taste symmetry, cut-off effects are still significant, and comparing to HRG is more subtle and 
requires further study, Fig. ||| In the left panels of Fig. || and || temperature is set by Sommer scale 
(n in this case). 

In part cut-off effects in the observables directly related to hadrons can be accounted for if 
the temperature scale is set with a hadronic observable. Following the insight by the Budapest- 



Wuppertal collaboration [12, 13] we use the kaon decay constant fx for this purpose. For the 
strange and light quark number susceptibility the results with fx scale are shown in the right panels 
of Fig. and |. 
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Figure 4: The strange quark number susceptibility for the HISQ/tree action at mi = m s /2Q with r\ (left) and 
fx (right) temperature scale. The solid curve in the right panel is from the HRG model. 
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Figure 5: The light quark number susceptibility for the HISQ/tree action at m/ = m s /20 with r\ (left) and 
fK (right) temperature scale. 



As one can see, for the strange quark number susceptibility the f% scale eliminates virtually 
all cut-off dependence and the data from different lattices are hardly distinguishable. For the light 
quark number susceptibility some residual cut-off dependence remains. This may be expected for, 
at least, two reasons: a) %l i s verv sensitive to taste symmetry breaking (even at our finest, N t = 
12 lattice the root-mean-squared pion mass is still about 200 MeV for the HISQ/tree action 
b) there is no a priori reason for the lattice spacing dependence in %u Xs to be the same as in fy. 

Next, we consider correlations of the electric charge and strangeness. With the r\ scale we 
observe substantial cut-off dependence, which is largely reduced with fx scale, compare left and 
right panel in Fig. ^[ At high temperatures these correlations are close to the value expected in 
the non-interacting ideal gas (Stefan-Boltzmann limit), while at low temperatures they are well 
described by HRG. 

Correlations of different quark numbers are also a convenient way to study the deconfinement 
aspects of the transition. At high temperatures such correlations are expected to be very small, 
while at low temperatures hadronic degrees of freedom give naturally rise to such correlations. In 
Fig. [7] we present our results for the light-strange and light-light quark number correlations. We 
divide them, correspondingly, by the strange and light quark number susceptibility in attempt to 
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Figure 6: Electric charge-strangeness correlations with r\ (left) and fx (right) temperature scale for 
HISQ/tree. The solid curve in the right panel is from the HRG model. 




Figure 7: u— and 5— quark and u— and d— quark 
correlations rescaled by the corresponding fluctua- 
tions with r\ scale for HISQ/tree. 



Figure 8: The trace anomaly for the asqtad 
and HISQ/tree action, m/ = m s /20. "s95p-vl" 
parametrization is explained in the text. 



get rid of the trivial mass effects. For the strange-light quark number correlations we see a good 
agreement with the HRG model, while the discrepancies with HRG predictions for Xud/Xi are 
likely due to cut-off effects. 

Comparing Figs. f|-|7] we can conclude that quantities, where pions do not contribute, can be 
reproduced quite well on the lattice, while those, for which pions give dominant contribution, still 
demonstrate noticeable cut-off effects, that are not completely removed by switching to a hadronic 
(fx) scale (at least, on the lattices typically in use, /V T = 6 — 12). 

4. Update on the trace anomaly 

The deconfinement transition is also often defined as the rapid rise in the energy density as- 
sociated with liberation of new degrees of freedom. The energy density and other thermodynamic 
quantities are usually obtained by integrating the trace anomaly (e — 3p)/T 4 . This quantity is under 
extensive investigation on the lattice. We present an update of the results reported in [jl|] in Fig. ||. 
The solid line is a parametrization of (e — 3p)/T derived in Ref. [[14]] that combines HRG result 
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at low T ( T < 160 MeV) with the lattice data 0] at high T (T > 220 MeV ). From the figure we 
conclude that the lattice results overshoot the HRG prediction for T > 160MeV since any cut-off 
effects could only decrease (e — 3p) /T 4 at low T. 

5. Conclusions 

In this contribution we studied the chiral transition in 2+1 flavor QCD close to the physical 
point. Using the 0(N) scaling we determined the chiral transition temperature in the continuum 
limit at the physical light quark mass to be 154(9) MeV. We also studied the deconfinement phe- 
nomenon in terms of the renormalized Polyakov loop as well as in terms of fluctuations and corre- 
lations of conserved charges. We concluded that it is difficult to study the deconfinement aspects 
of the transition in terms of the Polyakov loop, while correlations and fluctuations of conserved 
charges are much better suited for this purpose. If fx is used to set the scale we see a good agree- 
ment with the HRG model predictions for fluctuations and correlations that do not involve pions in 
their HRG expansion. For the quantities like %i and %ud, where the pion contribution is significant, 
cut-off effects are still too large for a meaningful comparison with the HRG model. 
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